Cosmological Imprints of a Generalized Chaplygin Gas Model for the Early Universe 
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We propose a phenomenological model for the early universe where there is a smooth transition 
between an early "quintessence" phase and a radiation dominated era. The matter content is 
modelled by an appropriately modified Chaplygin gas for the early universe. We constrain the 
model observationally by mapping the primordial power spectrum of the scalar perturbations to the 
latest data of WMAP7. We compute as well the spectrum of the primordial gravitational waves as 
would be measured today. We show that the high frequencies region of the spectrum depends on 
the free parameter of the model and most importantly this region of the spectrum can be within 
the reach of future gravitational waves detectors. 
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I. INTRODUCTION 

While our knowledge about the universe has improved 
over the last decades with the advent of new observa- 
tional data, there are several dark sides of the universe 
that have not been so far described from a fundamental 
point of view: what caused the initial inflationary era of 
the universe? what is the origin of dark matter? what is 
the fundamental cause of the current acceleration of the 
universe? Even though nowadays none of the previous is- 
sues have a satisfactory answer, a parallel approach, that 
can shed some light on the dark sides of the universe, is 
a phenomenological one or a model building strategy. A 
pod example in this regard is the inflationary paradigm 
l], where an inflaton field (or several scalar fields) in- 
duce the initial acceleration of the universe. Such a field 
would leave tracks on the evolution of the universe for 
example through the primordial power spectrum of the 
scalar perturbations 043] , extremely useful to constrain 
the inflaton field. It is precisely such an approach that 
we will follow in this paper. 

The main goal of this paper is to obtain a phenomeno- 
logically consistent model for the early universe (infla- 
tionary and radiation dominated epochs) by properly 
modifying the generalised Chaplygin gas (GCG) |g. A 
first attempt in this direction has been recently carried 
out in Q (see also 0) where a new scenario for the 
early universe was proposed. Such a scenario provides a 
smooth transition between an early de Sitter-like phase 



'Electronic address: mariam.bouhmadi@ist.utl.pt 
^Electronic address: chen@slac.stanford.edu 
tElectronic address: f97222009@ntu.edu.tw 



and a subsequent radiation dominated era. In that model 
Q , the matter content was given by a type of generalized 
Chaplygin gas for the early universe, with an underlying 
scalar field description. Here, we give a more realistic 
model where the early inflationary phase of the universe 
is more general than a de Sitter-like universe. We will 
show how the spectrum of the present model is not as 
red as the one presented in [|| and therefore more con- 
sistent with observations. 



More precisely, rather than having a de Sitter- 
like phase in the past, we will consider an early 
"quintessence" inflationary phase. This phase will be 
connected to a radiation dominated phase at later time. 
The model can be described through a scalar field or a 
Chaplygin gas inspired model. We will then analyze the 
possible imprints of such a gas in the primordial power 
spectrum of scalar perturbations and the power spectrum 
of the stochastic background of gravitational waves. 



The paper is organized as follows. In section II, we will 
present the model which is based on a properly modified 
Chaplygin gas that interpolates between a quintessence 
inflationary era and a radiation dominated period. In 
section III, we will show how this kind of gas can be 
modeled by a minimally coupled scalar field. In section 
IV, we constrain our model observationally and obtain 
the full spectrum of scalar perturbations numerically. In 
section V, we obtain the spectrum of gravitational waves 
using the method of Bogoliubov coefficient. Finally, in 
section VI we summarize our results and conclude. 
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II. THE MODEL BUILDING 

We generalize the model presented in Q by considering 
an inflationary period corresponding to a "quintessence" 
like behavior (described by a power law expansion) and 
followed by a radiation dominated epoch. The matter 
content of the universe can then be modeled a la Chap- 
lygin gas as 



P 



B 



A 



jl+/3 4(l+a) 



l/(l+a) 



(2.1) 



where A, B, a, are constants being A, B positive. 

The matter content is not interacting with any other 
fluid and therefore its energy density is conserved: 



P + 3H(p + p) = 0, 



(2.2) 



where the dot stands for derivative with respect to the 
cosmic time and p is the pressure of the fluid. If we 
substitute the energy density (|2.1[) in the conservation 
equation (|2.2p . we obtain 
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The above equation of state can be rewritten as 



V=\ P 



A 



1+0 -4(1 + a) 
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(2.3) 



(2.4) 



This equation shows clearly that we recover the model 
discussed by one of us in [(| for — > — 1 . We would like to 
highlight that the equations of state (|2.3j) - (|2.4|) have been 
previously analyzed in [|[ to study a possible interplay 
between dark matter and dark energy and therefore in a 
scenario completely different to the one we study here 1 . 

As the inflationary period takes place before the radi- 
ation dominated one, we deduce that the inequality 



B 



,4(l+a) 



A 



must hold at early time; i.e. at small scale factors. This 
will be the case as long as 

< 1 + -4(1 + a). (2.6) 

The inequality (|2.5j) implies a < ai, where 



a\ = 



l/[l+/3-4(l+a)] 



(2.7) 



The scale factor a\ represents the "border line" between 
two different regimes. More precisely 



^i/(i+«) 

a (l+P)/(l+a) 
BV(l+a) 



a <C a h 



a 3> a\. 



(2.8) 
(2.9) 



1 We thank Luis Chimento for pointing out this to us. 



There are two more conditions that should be imposed 
on a and to have a model that interpolates between 
an early inflationary phase of the type of quintessence 
and a radiation dominated phase at later time: (i) the 
energy density (|2 .8[) must induce a period of inflation 
and (ii) such a period of inflation should not induce a 
super-inflationary expansion; i.e. < H, and therefore a 
super-accelerating phase of the universe where the energy 
density (|2.8[) corresponds to phantom matter, with an 
energy density that grows as the universe expands. By 
combining these two ansatz with the inequality (|2.6[) . we 
can easily deduce that the set of allowed values of a and 
satisfy 2 



1 + a < 0, 
1 + /3 < 0, 
2(1 + a) < 1 



■13. 



(2.10) 
or (El 



A perfect fluid with an equation of state 
describe an early inflationary period as long as the con- 
ditions (|2.10[) arc satisfied. The universe will exit the 
inflationary epoch when the scale factor reaches the value 



a* 



2(1 + a) - (1+0) 
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(2.11) 

The last condition can be obtained by imposing the con- 
dition p + 3p = 0. Again this is in agreement with the 
results of for -> -1. 

Before concluding this section, let us analyze more 
closely what sort of inflation does the universe undergo 
in its initial stages. The Fricdmann equation with the 
matter content (|2.1[) is too difficult to be integrated an- 
alytically so we use an approximation where the energy 
density can be written as in Eq. (|2.8[) . Then the scale 
factor can be approximated as 



(2.5) a(t) 



1 + 
2(1 + a)' 



, A l/[2(l+ a )] 




2(l+a)/(l+/3) 



(2.12) 



where t is the cosmic time and k 2 — 8irG where G is the 
gravitational constant. As we can see from the previous 
equation the inflationary era of the universe is given by 
a power law expansion. 

In sections IV and V, we will see that the entire evolu- 
tion of the universe can be described by the appropriately 
modified Chaplygin gas (|2.4[) and a subsequent ACDM 
expansion, that is 



.1/(1+") 



l+i + o) . 



early-time 



(2.13) 



PMjr + Pmoijr + PA, late-time 



2 Having a negative value a is not that surprising (see for example 
Q). Indeed, it has been shown in ll(J| that type la Supernovae 
allow for negative a under the standard generalized Chaplygin 
gas. 
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FIG. 1: The blue curve corresponds to the scalar field, <j>, 
against x, where x corresponds to a power of the scale factor 
a, x = {B/A)a q , with q = 1 + j3 - 4(1 + a). The values 
xo and a;* correspond to the moments when the pivot scale 
ko — 0.002 Mpc -1 exists the horizon and the inflation ends, 
respectively. 
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FIG. 2: The blue curve corresponds to the scalar field po- 
tential V{(j>) against (f> (see Eqs. (|3,2[) - (|3,3[) ). where Vo = 
pjy B )-(i+/3)/[9U+«)]_ The values fa and fa corre _ 

spond to the moments when the pivot scale ko = 0.002 Mpc -1 
exists the horizon and the end of inflation, respectively. 



where ao is the current scale factor, p r Q, p m o an d Pa 
are the current energy densities corresponding to radi- 
ation, matter (cold dark matter and baryonic matter) 
and dark energy (modelled through a cosmological con- 
stant), respectively. On the other hand, at the radiation 
dominated epoch the energy densities (|2.13|) arc equal, 
consequently 



B = {p r0 a ) 



(2.14) 



The parameter A is related to the scale of inflation (see 
Eqs. HUH) and ([57f])). 



scale factor (see Fig. [T]): 
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(3.3) 



III. AN UNDERLYING SCALAR FIELD MODEL 



The inflationary dynamics of the model presented in 
the previous section can be described through a min- 
imally coupled scalar field, <fi, with a potential, V(<p), 
whose energy density and pressure read 



Zi CI Zi CI 

In the previous equations the prime stands for derivative 
with respect to the conformal time. This scalar field can 
be mapped to the perfect fluid with equations of state 
(|2.3[) or (|2.4j) ; i.e. = p and = p. Through this 
mapping we can write down <f> and V as a function of the 



where x = {B/A)a q and q = 1 + /3 — 4(1 + a). The scalar 
field starts with a negative value and it rolls down the 
potential as the universe inflates (cf. Fig.[2|). 

We can also consider the behavior of the equation of 
state parameter w, w = p^/p^, for the scalar field as a 
function of the scale factor (see Fig. [3]). At early times, 
the kinetic energy is negligible with respect to the po- 
tential, and therefore w — > — 1. Afterwards, the param- 
eter w increase gradually with time. The universe stops 
inflating when w reaches the value —1/3 and becomes 
radiation dominant when w finally reaches 1/3. 

At very early times where the scale factor is very small, 
it can be shown that the inflationary expansion follows a 
power law. Indeed, the scalar field and the potential can 
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be approximated by 

0(a) * ' ^ 1 ^ 
V(a) ~ 

00 = 



K V 1 + Q 
1 + 



6(1 + a) 




ln(a), (3.4) 
Al/(1+a)a _ (1+w/(1+a))(35) 



(3.6) 



of the slow-roll parameters e(a) and 77(a), where 
2 1 1 d 2 V 



1 

2^? 



1 dV 
V~d4 



n 



K V 



(3.12) 



and therefore 



V{4>) 



1 - 



1+/9 



6(1 + a) 



A 1 /'^") exp 



\l 7-J— (</> ~ 0o) 
1 + a 



At much later times where the scale factor becomes 
large, the universe is radiation dominant. It can be easily 
proved that for this period: 



(3.8) 
(3.9) 



0(a) ~ 


2 

- ln(a) + c, 
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V(a) ~ 
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nq 


In f 






+ /? 






+ a 




Consequently, we obtain 

7(0) ~ ex p [- 2k (0 -00] 



(3.10) 



(3.11) 



Since we know the form of the potential V(<j>) for this 
modified GCG model, we are able to find the evolution 




FIG. 3: The equation of state w for the modified GCG model 
with equation of state (|2.3|) or (|2.4p . The dashed red line 
indicates the end of inflation. We set 01 = —1.06 as an example 
here. 



The slow-roll approximation is valid when the condi- 
tions e < 1, i| C 1 arc satisfied, which means the poten- 
tial energy dominates over the kinetic term during the 
inflation era. In Fig. [4] we show the behavior of these 
functions in terms of the scale factor a. 



IV. PRIMORDIAL POWER SPECTRUM 

Inflation not only solve some of the shortcomings 
present in the big bang theory but also generates den- 
sity perturbations that seeds the structure of the present 
' universe. Those density perturbations have been con- 
\ strained through observations of the cosmic microwave 
background (CMB). In this section, we will constrain 
the model introduced in the previous section by using the 
measurements of WMAP7 [ll| for the power spectrum of 
the comoving curvature perturbations, P s = 2.45 x 10 -9 , 
and its index, n s — 0.963, where 



1 



d\nP s (k) 
din A; 



(4.1) 



These measurements correspond to a pivot scale fco = 
0.002 Mpc" 1 [Tlj]. 

Once the parameters of the model have been con- 
strained we will obtain the full spectrum of the scalar 
perturbations. 



A. Imposing Observational Constraints 

At very early times, the model introduced in Sect. II 
induces a power law expansion as shown clearly in 



2.0; 
1.5: 

1.0; 

0.5; 
0.0 L- 




1000 10 8 10 13 10 18 10 23 10 28 
a 

FIG. 4: The black line corresponds to e and the blue one to 77. 
The dashed red line locates the time when inflation ended. We 
also notice that the slow-roll conditions are no longer satisfied 
when the field is close to the end of inflation (a — a*). We 
set a — —1.06 as an example here. 
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Eq. (|2~T2|) (see also Eq. ([3J]) and Fig. HJ). On the other 
hand, the slow-roll approximation is valid for an inflation- 
ary power law expansion as the one we are considering 
(see for example Figs. [3] and U|). hi this regime, the power 
spectrum (for the comoving curvature perturbations) and 
the spectral index can be expressed as (see for example 



p, 



H 2 
2tt4> 

n s ~ 1 — 6e + 2r], 



(4.2) 
(4.3) 



Primordial Power Spectrum: Methods and 
Results 



The primordial density perturbation can be directly 
mapped to the comoving curvature perturbation, the last 
remains constant on large scale after exiting the Hubble 
horizon, as long as the perturbations are adiabatic. This 
is precisely the case, as we have a unique degree of free- 
dom corresponding to the generalised Chaplygin gas. 

The comoving curvature perturbation, s, is determined 
by the fluctuations of the generalised Chaplygin gas. The 
corresponding power spectrum for the field is [HQ 



at the horizon exit; i.e. at k = aH. Throughout the 
paper a dot stands for a derivative with respect to the 
cosmic time. For a power law expansion, i.e., a oc t 1 
where r is the conformal time, the spectral index reduces 
to n s = — (1 + 2/1). For the modified Chaplygin gas we 
are analyzing Z" 1 = (1 + /5)/2(l + a) - 1. 

With the previous inputs, we can proceed as follows: 
(i) the parameter B is fixed by the current amount of 
radiation in the universe as stated in Eq. (j2.14[) . (ii) for a 
given parameter a, the parameter j3 is fully determined 
by the measurement of n s and (iii) the parameter A 
is fixed such that P s = 2.45 x 10 -9 at the pivot scale 
ko = 0.002 Mpc -1 . Before moving forwards, we would 
like to highlight that the mode ko = 0.002 Mpc -1 exits 
the horizon well inside the inflationary era where both 
the power law expansion and the slow-roll approxima- 
tion are valid for the model (|2.4p . We found out that the 
inflationary scale, V = ^W+a) ( J 4/S)-(i+^)/?( 1 +«) ^ 
1.2 x 10 16 GeV, is almost of constant. 

The tensor power spectrum to the scalar power spec- 
trum ratio in the slow-roll approximation can be ex- 
pressed through the slow-roll parameter e 



Pt(k) 
P s (k) 



16e. 



(4.4) 



For power-law inflation in the slow-roll approximation, 
the spectral index is n s = — (1 + 2/1), as we already men- 
tioned, while e = 1 + 1/7. Because this GCG model be- 
haves like a slow-roll power-law inflation when the pivot 
kg exits the horizon, we can use these approximations to 
evaluate r at that time. Using once more WMAP7 data 
(n s k 0.963 at the scale fc = 0.002Mpc _1 ) we can find 
I. We obtain r = 0.296 at the pivot scale fc , which is in 
agreement with the WMAP7 constraints. We obtain the 
full spectrum of the gravitational waves on section V. 

At later times close to the end of inflation or when 
larger modes k exit the horizon, the slow-roll conditions 
are not fulfilled. This is a simple consequence of the be- 
havior of the GCG model that deviates from a power-law 
inflation at that time. Thus we cannot use the results in 
Eqs. (|4.2[) and (|4.3[) when inflation approaches a* where 
it halts. In this regime, we will rather calculate the spec- 
trum numerically. 



2ir 2 k- 3 P s (k) 



(4.5) 



where z = jj-. The field satisfies in the Fourier space 
the equation 0, [I[ 



k — — 

z 



vk = 0. 



(4.6) 



If the slow-roll conditions are satisfied, the variation of 
the field <fi an d the Hubble parameter H is much slower 
than that of the scale factor a. Therefore, the following 
approximation can be used: z" / z rj a" /a, at the lowest 
order of the slow-roll approximation. Consequently, the 
solutions of (14.61) read in this case H 



-ikr 



i 



i 

kr 



(4.7) 



where the Bunch-Davies vacuum has been imposed for 
modes well inside the horizon (aH <C k). We remind that 
the parameter r stands for the conformal time. Now, if 
we consider the super-horizon scale (|fcr| 1) and sub- 
stitute the approximation a ~ — 1/ (H r) into Eq. (|4.7|) , 
the solution becomes 



1 aH 



2k k 



(4.8) 



The corresponding power spectrum in this slow-roll ap- 
proximation is given in Eq. (|4.2p . These results remain 
valid on the next order of the slow-roll as shown for ex- 
ample in 0, 0] and therefore can be applied for an infla- 
tionary power law expansion as it is the case in the model 
we are analyzing at very early time. 

In order to obtain the full power spectrum P s (k), we 
use numerical methods instead of applying Eq. (|4.2I) for 
the reasons stated in the previous subsection. Therefore, 
we first solve the differential equation (|4.6p , and find the 
solution ^(ocross) calculated at the time when the mode 
k exits the horizon; i.e. k = a cmss H. Then through 
Eq. (|4.5j) we obtain the power spectrum P s {k)\k=aH ■ 

We proceed as follows. We separate Eq. f)4.6[) into two 
first order differential equations: 



X' = Y, 

Y' = -(k 2 -^)x, 



(4.9) 
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where we set X = v^. In order to solve the previous set 
of differential equations, we take the following actions: 

• It is easier to make a change of variable from the 
conformal time to the scale factor. Notice that z = 
a(j)/H can be fully determined in terms of the scale 
factor as <fi = aHd(f>/da and the scalar field <fi is 
fully determined as a function of the scale factor 
(see Eq. (pT2]>). 

• In addition, we have simply to use the set of values 
for the parameters of our model that have been 
deduced by imposing observational constraints on 
the model as discussed in the previous subsection. 

• Last but not least, we need to impose a set of 
boundary conditions. When the wavelength of a 
given mode k is much smaller than the Hubble ra- 
dius 3 k ^> aH, the effect of curvature can be ne- 
glected. Therefore, the result reduces to that of a 
flat Minkowski spacetime (when k ^ aH). So, the 
initial condition is 




2xl0- n i 



1 10 4 10 8 10 12 10 16 10 20 
k(Mpc _1 ) 



FIG. 5: Primordial power spectrum P s (k)\k= a H against k for 
six different values of a. The dashed black line is the pure 
power law inflation, and the vertical dashed red line locates 
the pivot k = 0.002Mpc"\ We can see that all these lines 
merge when small k; i.e. large scale, exits the horizon. The 
grey, violet, red, orange, green and blue curve correspond 
respectively to a = -1.1, -1.09, -1.08, -1.07, -1.06, -1.05. 




— iki 



for fcrl > 1, 



(4.10) 



which corresponds to the Minkowski vacuum at 
very early time. We change the variable from the 
conformal time r to the scale factor a through the 
relation 4 



a(r) 



1 + /3 
2(1 + a) 



- 1 



8ttG 



^l/[2(l+a:)]. 



2(1 + Q) " 



(4.11) 



Following a similar procedure, we can see if the spectral 
index, n s (defined in Eq. (|4.1[) h depends on the scale. 
Our results are shown in Figs. [5] and [5] 

For lower modes, i.e., those that entered the hori- 
zon very recently, the power spectrum is independent 
of the parameter 5 a. It is only for modes satisfying 
10 5 Mpc _1 < k that we can start to see some dependence 
of P s on a. As should be expected, the power spec- 
trum has a constant slope; i.e., a constant spectral index, 
for lower k when the modified Chaplygin gas induces a 
power-law expansion. Indeed, the results for power-law 
matches very well those of our model for these modes 
(cf. Figs. [5] and [6]). However, the spectral index of the 
primordial power spectrum is not scale independent as 
Fig. [5] clearly shows. The deviation of n s from the scale 



3 This condition is fulfilled by any mode (on the past) thanks to 
the inflationary mechanism. 

4 Notice at this regard that the expansion of the universe for modes 
with a wavelength much smaller than the Hubble radius is well 
approximated by a power law. 

5 The other parameters of the model are fixed for a given a as 
explained in the previous subsection. 



independence starts when the slow-roll condition ceases 
to be valid and therefore approximately a bite before the 
end of inflation. 

Before concluding this section, we recall that the cur- 
rent CMB observations cover roughly the range of scale 
k « 0.0002Mpc~ 1 to k w 0.2Mpc _1 In this observ- 
able region, it is hard to separate a modified GCG model 
for different values of a or even from a power-law infla- 
tion. Could currently running observations or future mis- 
sions help us in this regard? For example, it is expected 
that the maximum scale I in CMB observations from the 
Planck mission would reach £ max « 2000. This scale 
(■max could be translated into a k max scale. Even though 
the relation between i and k is not a one-to-one rcla- 




k (Mpc 



FIG. 6: The spectral index n 3 against k. The red dashed line 
locates the pivot scale ko = 0.002Mpc _1 . The spectral index 
n s is approximately a constant at early time. The grey, violet, 
red, orange, green and blue curve correspond respectively to 
a = -1.1, -1.09, -1.08, -1.07, -1.06, -1.05. 
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tionship, we can still find the main contribution for each 
mode: k (r — r*) ~ £ [lH, where r — is the elapsed 
conformal time since the last scattering surface until the 
present. Here we use the concordance ACDM model to 
calculate (to — r*) = J°° da/a 2 H numerically. We find 
that the maximum scale corresponding to l ma x ~ 2000 
is k m ax ~ 0.138Mpc _1 . Consequently and unfortunately, 
the answer to the question raised in this paragraph is neg- 
ative, as is clear from Figs. [5] and O We will next show 
that by looking at the high frequency spectrum of the 
gravitational waves, it is possible to separate a modified 
GCG model for different values of a. Those frequencies 
might be within the reach of future grav itational- wave 
detectors such as BBO and DECIGO jlj ]. 




a 



GRAVITATIONAL WAVE SPECTRUM 



FIG. 7: The integration method we use for a given mode 
k. The dashed line corresponds to (aH) 2 , and the solid line 
corresponds to a" /a for a modified GCG model. 



The early inflationary era not only leaves imprints on 
the scalar cosmological perturbations but also creates 
a fossil of gravitational waves. In this section, we will 
analyze the possible imprints in the power spectrum of 
the stochastic background of gravitational waves, for the 
model presented in section IF In this case the background 
evolution of the universe, till nowadays, is described by 
the matter content in Eq. (|2.13[) . This analysis is quite 
important as it can shed some light on the inflationary 
scenario behind the early accelerating phase of the uni- 
verse and its transition to the subsequent radiation era. 



A. Method 



(X - F) 2 /4, where X, Y fulfil 
( X' = -ikY 

\Y' = -i (k 2 -%)X- (5 - 2) 

The coefficient gives the number of gravitons, Nk, for 
each mode k created during the evolution of the universe, 
where N k (r) = \(3 k (r)\ 2 . 

To calculate the GWs spectrum (|5.1|) , we have to solve 
the differential equations (|5.2[) numerically and use ap- 
propriate initial conditions for X{ti) and Y{ri). We pro- 
ceed as follows: 



In order to obtain the spectrum of the gravitational 
waves (GWs), we will use the method of Bogoliubov co- 
efficients. To our knowledge, this method was first devel- 
oped in [l5rll3] and later applied in fl8i - l2lj |. Bogoliubov 
coefficients describe how the vacuum changes as the uni- 
verse expands. In particular, one of these coefficients, 
which we will denote (3k, gives the number of gravitons 
created. Furthermore, it can be shown that the dimen- 
sionless relative logarithmic energy spectrum of t he g rav 
itational waves, Slew, at the present time reads 20. [21 



• The expression for a" /a can be obtained from the 
Friedmann equation and the conservation law: 



a" k 2 2/ 
a o 



(5.3) 



and therefore can be written in terms of the scale 
factor, a, in this modified GCG model as follows, 



n GW (uj,To) 



1 dp GW 
Pc{t ) dlnw 



hn 2 



(To) = 3^H 2 (r f 4 ^ To) - 
(5.1) 

The parameter pew is the energy density of GWs and 
u the respective angular frequency; p c and H are the 
critical density of the universe and Hubble parameter, 
respectively, evaluated at the present time. 

In summary, the present value of the Bogoliubov coef- 
ficient, /3k, determines the power spectrum of the grav- 
itational waves. This parameter can be determined in 
terms of two continuous functions X, Y such that \f3 k \ 2 = 



a 



-a 

A 
1 2 

-a 



111 

l+c* 
B 



A 



X 

a/(l+a) 



Pmo(f) +4p A 



, early time (5-4) 
late time 



• At very early time, the model we are analyzing be- 
haves like a pure power-law inflation. In that case, 
the set of differential equations (|5.2[) have an an- 
alytical solutions [2l|, which we will use as initial 
conditions for X(r) and Y(t) in our numerical in- 
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tcgration 6 : 



Y( n ) = -i 



— kTilT 



(5.5) 
(5.6) 



Before proceeding, let us highlight that the general 
solution of Eq. (|5.2j) for a power-law inflation in- 
cludes the first and the second kind of Hankel func- 
tions [22[ . But here we adopt just Hankel function 
of the first kind because the solution should reduce 
to the Minkowski spacetime solution (|4. 10[) at very 
early time when aH <C k Q ■ 

• Similarly to the scalar perturbations, it is simpler 
to solve the set of equations (|5.2[) in terms of the 
scale factor. With this in mind and in order to 
apply the boundary conditions (|5.5[) and (|5.6[) , we 
use the relation between the conformal time and 
the scale factor stated in Eq. (|4.11j) . 

• An important issue is the range of integration of 
Eqs. (|5.2[) . In principle, the integration must be 
done from a very early time when the mode is well 
inside the horizon on the inflationary era to the 
present time. However, there is a way of reducing 
the computing time without affecting in practice 
the results: for a given mode fc, we start integrat- 



ing from 0.01a crosSl until 100a c 



where 



and i 



correspond to the scale factor where the 



mode k exits the horizon and reenters the horizon, 
respectively. Our numerical integration will be car- 
ried out for modes k ranging from the maximum of 
the potential a" /a until the minimum one. The 
last one corresponding roughly to the modes that 
are reentering currently the horizon. The methods 
is schematically shown in Fig. [7] 

As can be noticed from Fig. [3 there is an intersec- 
tion point in the potential a" /a which takes place 
roughly at 10 50 . Given that we have set the cur- 
rent scale factor to ao = 10 58 , the intersection point 
corresponds roughly to the time of Big Bang Nucle- 
osynthesis (BBN). On our model, that point corre- 
sponds to the transition between the modified gen- 
eralized Chaplygin gas and the ACDM model. For 
modes that enters the horizon after the BBN, we 
will have to integrate Eqs. (|5.2p twice as we will be 
dealing with two different potentials as shown in 
Eq. (|5.4[) . For the second integration, we use the 
final values obtained in the first integration as the 
initial values for the functions X and Y . 




Notice that the conformal time is negative. 



1Q 36 1Q 45 



FIG. 8: Evolution of the graviton numbers, \/3k\ 2 = N k . The 
integration has been carried out for six different values of 
frequencies. We can see that for each mode k, the gravitons 
have been created mostly during the inflationary era. For this 
plot we fixed a = —1.06. 



The integration will be performed assuming: f2 r0 = 
8 x ICT 5 , fi m0 = 0.24, a = 10 58 and H = 71.0 
km/s/Mpc. In addition, the parameters of the 
model are fixed as explained in section II. 



B. Results 

In this subsection we present our results of the power 
spectrum of the gravitational waves for the model pre- 
sented in section II. 

As we explained before, the Bogoliubov coefficient, 
/3fc, gives the number of graviton created in each period, 
\Pk\ 2 — Nk(r), and therefore we can calculate the time 
evolution of during the expansion of the Universe in 
this GCG model. In Fig. [5] we present an example of our 
results where we set a = —1.06. The plot corresponds 
to six different wave numbers k: k = 1.9 x 10 _3 Mpc _1 , 
k = 1.9Mpc _1 , k = 1.9xl0 3 Mpc~\ k = 1.9xl0 6 Mpc _1 , 
k = 1.9 x 10 9 Mpc~\ and k = 1.9 x 10 12 Mpc _1 , or dif- 
ferent frequencies lo = k/cto. 

Similarly we can obtain the spectrum of the gravita- 
tional waves as shown in Fig. [3] This figure is quite en- 
lightening. It shows that for large frequencies the larger 
| a | is, the larger is the dimensionless logarithmic energy 
spectrum of the GWs. On the other hand, the larger \a\ 
is, the more the spectrum is shifted towards larger fre- 
quencies. This is a simple consequence of the fact that 
a larger \a\ implies a larger maximum of the potential 
a" /a. In summary, an increase in |a| implies two things: 
an upward shift and a rightward shift of the spectrum. It 
is also worthy of notice that the plateau of the spectrum 
merges at middle/low frequencies for different values of 
a. This is not surprising as the energy scale of inflation, 
Vq, is almost fixed and independent of a in our model, 
where Vq ~ 1.2 x 10 16 GcV (cf. the previous section) 
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where ho = 0.71. All these constraints are fulfilled by the 
power spectrum shown in Fig. [5] 

This model is a clear example of how the detection 
of GWs can be extremely helpful to distinguish among 
several inflation models and a complementary tool to the 
study of primordial spectrum of the scalar perturbations. 
Indeed, while our model is indistinguishable for any value 
of a at low k and ui, on the spectrum of scalar and tenso- 
rial perturbations, respectively, it should be distinguish- 
able at high k and lo. Most importantly, the high fre- 
quency regime is within the reach of future grav itational 
wave detectors such as BBO and DECIGO 1141 1 . 



FIG. 9: The gravitational wave spectrum Q.gw against the 
frequency to for different value of a in this GCG model: the 
blue line refers to a — —1.06, the red one refers to a = —1.05, 
and the grey one to a = —1.04. 

and it is precisely the value of Vb that shifts vertically 
such a plateau [2l| . The parameter Vb is fixed mainly 
by the measurements of P s and n s . The fact that n s is 
much more constrained by the latest WMAP7 data than 
the previous one, WMAP5, explains the changes in the 
plateau of our Fig. [9] and Fig. 4 of Ref. [2l[, where on 
the last mentioned figure different values of n s , or equiv- 
alently different values of the power law universal ex- 
pansion I (in terms of the conformal time) are assumed. 
Therefore, the spectrum of GWs at very low frequencies 
is insensitive to the value of a, because for larger fre- 
quencies the spectrum already merged for different a's 
whereas for smaller frequencies the universe is described 
by the concordance model ACDM (at time the modes 
reenters the horizon), which is independent of the pa- 
rameter a. 

There are some observational constraints on the upper 
limit of the energy spectrum £l G w 0] : 

• Constraint from CMB: 
hln GW (u} hor ,T ) < 7 x 10 -11 , 
ioTtUhor = 2x 10 _1 ho rad/s. 

• Constraint from timing observations of millisecond 
pulsars: 

t ) < 2.0 x 10~ 8 , 
for uj pu i = 2.5 x 10~ 8 rad/s. 

• Constraint from doppler tracking of the Cassini 
spacecraft: 

hln GW (uJcas,To) < 0.014, 
for LOcas = 7.5 x 10~ 6 rad/s. 

• Constraint from LIGO: 
h 2 Sl GW (u>,T ) < 3.4 x 10- 5 , 
for a few hundred rad/s. 

• Constraint from BBN: 

J^ 3 Q G w(w, T )duj/uj < 5.6 x 10~ 6 , 
where w„ w 10 -9 rad/s, 



VI. CONCLUSIONS 

We present a model that attempts to fuse the infla- 
tionary era and the subsequent radiation dominated era 
under a unified framework so as to provide a smooth tran- 
sition between the two. The model is based on a modifi- 
cation of the generalized Chaplygin gas. More precisely, 
it interpolates between a "quintessence-like" , or power 
law, expansion and a subsequent radiation dominated 
universe. Such a gas fulfills an equation of state (|2.3j) 
or (|2.4|) which has been previously used in the literature 
as a means to explain a possible interplay between the 
dark sectors of the Universe Q. Our attempt is there- 
fore in a different context. We have shown how such a 
model has an underlying scalar field description, where 
the scalar field starts rolling down a potential hill where 
the extreme slow-roll approximation is valid until a bite 
before the end of inflation. 

We have obtained the full power spectrum of the scalar 
perturbations and constrained the model using the latest 
WMAP7 data. From our analysis it turns out that the 
model is indistinguishable from a power-law expansion 
at low wave numbers k (cf. Figs. [5] and [6]) even if we 
vary the only degree of freedom present on the model, 
the parameter a. Notice in this regard that the rest of 
the parameters of the model are fixed by observations as 
explained on Sect. IV. 

We have completed our analysis by obtaining the spec- 
trum of the gravitational waves (cf. Fig. [SJ where we 
have used the method of Bogoliubov coefficients fl5l - [2~0 ] . 
It turns out that our model is within the reach of fu- 
ture gravitational-wave detectors like BBO and DECIGO 
Most importantly, for those frequencies it appears 
possible to distinguish this modified GCG model among 
different values of a, while this seems unattainable using 
the power spectrum of the scalar perturbation with the 
present or near future observations. 

Last but not least, we have shown as well how the 
scalar perturbation spectrum of the present model is not 
as red as the one presented in Q for an alternative GCG 
model for the early universe. In comparison, the model 
presented here is more consistent with observations. 
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